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Abstract 

We investigate the spin-orbit induced spin-interference pattern of ballistic electrons travelling 
along any regular polygon. It is found that the spin-interference depends strongly on the Rashba 
and Dresselhaus spin-orbit constants as well as on the sidelength and alignment of the polygon. We 
derive the analytical formulae for the limiting cases of either zero Dresselhaus or zero Rashba spin- 
orbit coupling, including the result obtained for a circle. We calculate the nonzero Dresselhaus and 
Rashba case numerically for the square, triangle, hexagon, and circle and discuss the observability 
of the spin-interference which can potentially be used to measure the Rashba and Dresselhaus 
coefficients. 
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I. INTRODUCTION 



In recent years, much research has been devoted to study means to manipulate the spin 
of electrons. A prime candidate for spin-manipulation is the Rashba spin-orbit interaction 
which arises in certain semiconductor quantum wells such as InGaAs heterostructures 
as a result of structural inversion asymmetry. The strength of the spin-orbit interaction de- 
pends on the design of the heterostructures jjj and can be dynamically altered by applying 
a gate voltage 0] , hence offering a way to manipulate the spin in a controlled fashion. 
Another spin-orbit interaction present in for instance GaAs is the Dresselhaus spin-orbit 
interaction [5] which stems from the absence of bulk inversion symmetry. Recently j(| it has 
been shown for bulk GaAs structures that induced uniform strain affects the Dresselhaus 
spin-splitting hence giving another way to engineer spin manipulation, 
n this article, which is a justification and generalization of our recently introduced proposal 
3], we will further develop our idea to study spin manipulation by means of the Rashba and 
Dresselhaus interaction via spin-interference effects on the conductance. We will treat the 
Rashba and Dresselhaus spin-splitting on an equal footing and we will derive the interference 
pattern resulting from their combined effect on the spin of the electron. 
For ring-structures the interference effect due to spin-orbit interaction has been well inves- 
tigated 8, and for square loops it has recently be shown |13( that interference 



due to Rashba spin-orbit coupling can lead to electron localization. We will discuss the 
interference pattern due to the Rashba and Dresselhaus spin-orbit interaction of electrons 
travelling along polygons. In this way we are able to propose a scheme how to experi- 
mentally obtain the Rashba and Dresselhaus parameters from the spin-interference. The 
electrons will be treated as ballistic wave packets with a spin travelling both clockwise and 
counter-clockwise along the sides of a regular polygon which is positioned in the x-y-pl&ne, 
the 2DEG of a semiconductor heterostructure. Since the spin-orbit interaction depends on 
the travel direction, the spins of the electrons will rotate along different directions for each 
side of the polygon with a precession angle that depends both on the sidelength as well as 
on the strength of the Rashba and Dresselhaus spin-splitting. For certain precession angles 
the counter-clockwise and clockwise spin waves interfere completely destructive or construc- 
tive after a full rotation, which corresponds respectively to a large and small conductivity. 
We include an Aharonov-Bohm (AB) phase difference between the counter-clockwise and 
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clockwise spin waves and derive that for an unpolarized beam the spin-interference affects 
the amplitude of the AB oscillation. The procedure we use to calculate the spin-interference 
transforms the spin by a series of quantum-mechanical rotations corresponding to the various 
sides of the polygon. 

II. INTERFERENCE INTENSITY 

Consider the n-sided regular polygon of Fig. ^ which the electron traverses both clockwise 
and counter-clockwise simultaneously. Let us concentrate first on the counter-clockwise 



B 

FIG. 1: Polygon with n vertices. The effective magnetic field direction changes as a function of 
the side k. 

travelling part. The spin-orbit interaction manifests itself as an effective magnetic field 
about which the spin precesses. This effective magnetic field is constant along a side but 
differs for different sides. We assume for the moment that the electron will start at the 
beginning of the k = side indicated by the in Fig. ^ which we will call the origin. We will 
denote the quantum mechanical rotation of the spin along a side k with Rk and the precession 
angle (of which Rk is a function) with s. Suppose that the initial spin polarisation direction 
is along the z-axis. We thus have a J 2 and J z eigenstate |§,77i) = \m). The wavefunction 
after travelling along the side k then reads 



where D^^Rk) is the rotation matrix |l4j corresponding to the rotation operator Rk. Using 
the group multiplication law for the rotation matrices repeatedly, we find for the spin wave 




y 




(1) 



m' 
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function after completing the polygon counter-clockwise 

A\m)=J2D 1 J l UA)\m>), (2) 

m' 

where the rotation A is given by 

A = R n -i • ■ • R 
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The rotation A is characterised by the rotation direction unit vector £ and angle n( which 
we will calculate analytically later for some limiting cases. The clockwise going spin wave 
function reads 

A" 1 \m) = Y.D 1 ^' 1 ) W) = E^w(^) K> » (4) 

m' m' 

where the last step follows from the unitarity of the Wigner D-functions and * denotes the 
complex conjugate. The total final state becomes 

\^) = ^A\m) + ^=A- 1 \m), (5) 

yielding for the interference intensity in the origin 

\m\ 2 = 1 + Re (DU 2 m (A*)) . (6) 

We now generalize by starting with a spin polarised in an arbitrary direction, which we 
however can decompose into a linear combination of spin-up and spin-down along the original 
z-axis 

|/m")=a|T) + &|l>- (7) 

Furthermore, we include the effect of a magnetic field causing an Aharonov-Bohm phase 
difference 5 between the counter-clockwise and clockwise travelling wave (we ignore the 
effect of the magnetic field on the spin). We find for the total final state 

m = e lS ^A |/ m") + ±=A~ X \/ m") , (8) 



and for the interference intensity 

|^(0)| 2 = H 2 (l + Re (e^Dl^A 2 
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+ \b\ 2 il + Re(e t5 DK,(A 2 )jj (9) 
+ 2Re (a*b (e~ iS D^* Ll (A 2 ) + e i5 Df*_ L (A 2 )^ . 



In the following we will frequently make use of Euler parameters to express rotations, for a 
discussion of Euler parameters see, for instance, The interference intensity in terms of 
Euler parameters reads 



|^(0)| 2 = 1 + e cos<5+ (H 2 - |6| 2 )e 3 sin5 
+ 4 [Re(a*b)e 1 + Im(a*b)e 2 } sin 5. 

If we average over the spin polarisation the intensity reduces to 



(10) 



|^(0) | 2 = l + e (A 2 )cos6. 



(11) 



In the following Section we will discuss the various contributions to the rotation A 2 and seek 
a closed expression for it. 

III. ROTATION 

The important spin-orbit contributions to the Hamiltonian are the Rashba and linear 
Dresselhaus terms, for a [001] 2DEG given by, respectively, 



H R = a R {xk y - yk x ) 
H D i = ac m (-xk x + yk y ) 



(12) 



Their effect on the spin can be described as a rotation about the sum of the effective magnetic 
fields of the individual terms, the (dimensionless) direction of which is given by the vector 



( sin/3 \ 



B = s R 



\ 



cos (3 




( - cos p \ 



+ s m 



\ 



sin (3 




(13) 



/ 



which follows from Eq. ([12)1 by expressing the wavevector in the angle (3 with the x-axis. 
The parameter Sj, i = R, Dl is the angle through which the spin would rotate about the 
effective magnetic field direction Bi in the absence of the other. Hence Sj determines the 
relative weight of that particular term in the sum. The precession angle about the total 
effective magnetic field is the modulus of the vector Eq. (fT3*j) 



\B\. 



(14) 
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The Euler parameters then read 

e = cos — — 

B <15) 
e = B sin ^— L . 

2 

For a polygon, the angle (3 assumes only the discrete values (3 k = 2irk/n where n is the 
number of vertices of the polygon and k ranges from to n — 1. Notice that we restrict 
for the moment 1 edge of the polygon to be positioned along the x-axis (as in Fig. [IJ, a 
condition on which we will comment extensively later on. In the expectation that rotations 
along successive polygon edges will partially cancel each other we decompose the rotation 
along a side into Euler angles (<f,9,ip): Rk = R z (<p)Ry(9)R z (i(i) (We use the socalled zyz- 
convention). We find for the Euler angles 

9 = sign(B)s((3 k ) 

¥ = -ip (16) 
s R sin (3 k - s D1 cos (5 k 
-s R cos Pk + s m sm(j k 

where we have included the overall sign of B in the Euler angle 9 since it is divided out in 
the expression for the other Euler angles. The Euler angles become particularly simple for 
sm = ±s R for then ip = —Lp = 7r/4,37r/4, 9 = v^s^sin (3k =F cos (3k), respectively, and the 
total rotation Eq. © becomes 

A = R z Ry (^V2s R J2sm(3 k T cos/3 fc j R z = 1. (17) 

Hence we always have a maximum for % = ±sr. 



A. am = or olr = 

Consider now the limit odi = 0. We find for the Euler angles 

{<p,e,i/>) = (j3 k ,s,-/3 k ), (18) 

where we defined a new variable 

% - — , (19) 

n 



and where we have written s for sr. The rotation Rk we may then write as 



Rk = Rz(Pk)Ry(s)R z (—f3k)- 



(20) 



The product of two adjacent rotations in Eq. Q takes the form 



Rk+iRk = Rz{l3k+i)Ry{s)R z {—l3i)R y {s)R z {—(3k 



(21) 



The total rotation Eq. (jSJ) then simplifies to 

A = R n ^ ■■■R Q = R z ((3 n ) [R z (-(3 1 )R y (s)f 

where we introduced a new rotation 



(22) 



R^C)=R z (-p 1 )R y (s). 



(23) 



The interpretation of the new rotation Eq. ()23j) is straightforward: first a rotation about y 
and then, since the magnetic field rotates counter-clockwise, the spin vector rotates clockwise 
about z. In terms of Eulerparameters the rotation Eq. ()23|) reads 

(24a) 
(24b) 
(24c) 
(24d) 



Pi 8 C 

e = cos — cos - = cos - 
u 2 2 2 

. Pi . s , . C 
ei = sin — sin - = f x sin - 
2 2 2 

Pi ■ s , . ( 
e 2 = cos — sin - = f v sin - 
2 2 v 2 



. C 



sin — cos - = iz sin -, 



yielding for the rotation-direction 



1 - cos 2 y cos 2 | 



V 



sin — sin | » 

8i • s 

cos <y sin 2 
sin y cos | j 



(25) 



The Euler parameters of the square of the total rotation Eq. (JHJ) (i.e. A 2 ) are 



e = cosn£, 
e = £sinn£. 



(26) 
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From Eq. (jllj) we can directly read off the condition for a maximum respectively a minimum 
in the interference intensity in terms of the Euler parameters 



e 
e = 



= 1 
-1. 



(27a) 
(27b) 



Eqs. (|27a|) respectively ()27b|) together with Eq. specify the angles ( for which we 

have a maximum respectively a minimum. Combining Eq. (|27aj) with Eq. (J24a|) we find a 
maximum intensity if 



cos s 



COS ^ + 1 
n 

cos 2 ^ 



1, 771 £ 



From Eq. (j27b|) together with Eq. 



(2m+l)7r . -i 
COS h 1 



aj) we find a minimum intensity if 

71 — 1 



COS s 



cos 2 ^ 



— 1, m G 



(28) 



(29) 



The restrictions for m follow from coss G [—1, 1]. 

From Eq. ()24a)l we see that if (and only if ) s — > s + 2ir, then ( — > £ + 27r. The rotation 
A 2 is characterised by n2£ which transforms then as n2( — > n2( + n47r = n2(. Thus A 2 
is periodic in s with period 27r, ergo is the interference intensity. That is the reason why 
we could rewrite the conditions to get extrema for coss/2 into conditions for coss without 
loosing any information. 

Apart from the round trip interference which happens at the entrance lead of the polygon 
we could also imagine a lead halfway of the polygon causing interference in the forward 
direction. It can be readily observed that the forward interference intensity is given by the 
same expressions as the backward interference, except for the substitution A 2 — > A. The 
intensity halfway attains maxima and minima, respectively, for 

n~ 



(2m + l)lT 7T s 

cos = cos — cos-, mG 

n n 2 



2mTX 



1 4J 
n 

'4 



(30a) 

it s r ni 

cos = cos — cos-, mG 1,— . (30b) 

n n 2 L 4J 

Notice that for n odd we cannot rewrite the condition for coss/2 in a condition for coss 

since the period in s is 4n for n odd. 

If the starting point of the spin wave function is not the origin but, say, p radians removed 
from the origin (on the k = side), the total rotation becomes 



A' = R y (p)AR-\p), 
8 



(31) 



which is merely a similarity transformation and transforms the spin to a different basis. 
Since we average over the spin a similarity transformation does not change the result for 
the intensity. If the polygon has not one side aligned along the ^-direction but tilted at 
some angle x then we have the substitution (3k — ► (3k + X- Again, the total rotation will be 
changed by a similarity transformation 

A' = R z (x)AR~ 1 (x), (32) 

unaltering the results for the intensity. Things change however if both the Rashba and 
Dresselhaus interaction play a role, as shown in Section fill Bl since then adjacent rotations 
do not cancel each other neatly. 

The Rashba 5-field becomes the linear Dresselhaus 5-field if we make the substitution 
Pk ~ > — 7r/2. The total rotation Eq. (j3J) once again changes by a similarity transformation 



A' = R z (-n/2)AR;\-ir/2). 



(33) 



Hence the derived results hold as well in the limit or ; = 0, am ^ 0. 

We will now rewrite Eq. (fTTj) in terms of physical variables. According to the spin 
precession angle s in a one dimensional ballistic channel with Rashba- SO coupling is given 
by 

2am* 

s = _ x L, (34) 

with a the Rashba coefficient, m* the effective electron mass and L the sidelength of the 
polygon. Using Eq. ()24aj) to express ( in s we find for the averaged interference intensity 
Eq. (HQ) 



\m\ 2 

1 + cos I 2n arccos 



7T 

cos — COS 
n 



( 



2am* 1 



L 



cos 5, 



(35) 



V h 2 2 

where we substituted Eq. (134)1 for the precession angle s. An interesting limit is to let the 
number of vertices n go to infinity and at the same rate reducing the sidelength L, thereby 
obtaining a circle. From Fig. ^ we derive the relation 



7T 



sin 



n 



L 

2R' 



(36) 
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with R the circumradius. We can express n in terms of L by using Eq. (|36|). and taking the 
limit L — > we find for the intensity 

KKOlF Icircle= 

/ / 7w~v\ r ( 37 ) 

1 +cos I 2wd 1+ ( — ^— Rj cos 5. 

The intensity halfway is obtained from Eqs. ()35l 137)1 by dividing the argument of the cosine 
by 2. 



B. or, odi 7^ 

For general or, ckdi the above analysis does not lead to a closed expression for the total 
rotation. The general case can however be easily calculated numerically by multiplying rota- 
tion matrices and plugging the result into Eq. The result obtained for the interference 
intensity in the origin for the square is shown in Fig. El Fig. El shows the interference pattern 
for different alignments of the square in the xy-plane from to 7r/4 angle with the i-axis. 
The pattern is symmetric with respect to the offset angle ir/4 because of the symmetry be- 
tween the Rashba and linear Dresselhaus interaction and the pattern repeats itself of course 
after offset n/2. We observe as derived in Eq. (|T7)) the maxima lines or = ±o:di and the 
symmetry between Rashba and linear Dresselhaus. The maxima and minima along the lines 
sdi = 0, sr 7^ and sr = 0, sdi ^ follow from Eqs. (|28)l and (|2U)) with n — 4. They are 
given by s max = 0,tc,2tt & s min = arccos(l — \/2),27r — arccos(l — \/2). As we discuss in 
the following Section, the strong dependence on the precession angle and alignment can be 
exploited to observe the spin-interference. 



IV. DISCUSSION 



As a first order approximation the conductance is inversely proportional to the interfer- 
ence intensity in the origin |17|. hence any alteration of the spin-interference should manifest 
itself as a change in the conductance. Furthermore, as observed from Eq. pijl . the spin- 
interference affects the amplitude of the AB oscillation which is an experimentally measur- 
able quantity in solid state devices (la ]. In the following discussion we will therefore assume 
that the interference intensity can measured and we will focus on the spin-interference itself 
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and we will stipulate how to obtain the spin-orbit parameters from it. 
From Fig. |21 it is observed that the spin-interference varies as a function of the precession 
angle and the alignment of the polygon (a square in this case). The precession angle is 
a function of the sidelength of the polygon as well as a function of the Rashba and Dres- 
selhaus coefficients «j, i = R,D1, see Eq. (|34j) . With both the Rashba and Dresselhaus 
coefficients experimentally accessible parameters and the sidelength and alignment at choice, 
the spin-interference can be experimentally changed. In turn, the known dependence of the 
spin-interference on the sidelength and the alignment can be used to find the Rashba and 
Dresselhaus coefficients. As an illustration we have plotted in Fig. El the spin-interference 
of the square as a function of the sidelength and the alignment for realistic values for the 
Rashba and Dresselhaus constants. Varying the sidelength corresponds to moving along a 
line through zero with slope determined by the ratio of oir and am in one of the figures 
of Fig. El From measurements with different sidelength L as in Fig. El it is possible to 
draw contour lines of the measured interference intensity through which the lines should 
pass. Multiple measurements will single out a specific line hence giving the ratio of the 
coefficients. The order of magnitude of the coefficients can be quickly determined from the 
number of extrema within a certain range of the sidelength L. To find the coefficients indi- 
vidually, it will be necessary to fit the interference pattern with our model. In that respect 
it is promising that for different values of the coefficients the spin-interference displays very 
distinct behavior, which is illustrated in Fig. El which plots the interference pattern for two 
values of cxr. 

Until now we have constricted the discussion to square patterns but the type of polygon 
can be exploited as well for the determination of the spin-orbit parameters. In Fig. 0] we 
have plotted the interference pattern for the triangle (upper left), square (upper right), and 
hexagon (lower left), all aligned along the x-axis. The spin-orbit parameters can then be 
determined by measuring the interference intensity as a function of polygon which, by com- 
parison with our model, will yield the parameters. Especially the interference patterns for 
polygons with low number of vertices are distinctive. With increasing number of vertices, 
the polygonal interference pattern converges rapidly to that of a circle shown in the lower 
right of Fig. |U (which is actually a 1000 sided polygon), see below. 

In a real sample there might be two or more current leads attached to the polygon giving 
interference effects at all of them. However, the interference at the lead(s) different from 
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the one through which the electron enters will be smeared out by placing a large number of 
polygons in a row since the interference of the space-part of the wavefunction at the lead(s) 
is uncorrelated between polygons. This has been nicely shown in ^| to be the case for the 
AB h/e oscillations. 

In general, another spin-orbit contribution to the Hamiltonian, the cubic Dresselhaus term 
given by 



can become important for certain heterostructures. This term can be readily included nu- 
merically. We observe however that the cubic Dresselhaus term vanishes identically for the 
square aligned along the a:-direction and its effect can accordingly be minimized. Further- 
more, from simulations it is observed that, for reasonable values of am, hs effect on the 
spin-interference for the square aligned along the x — y-direction is minimal. Notice from 
Fig. El that for these two alignments of the square the spin-interference displays distinct 
features hence the alignment would still count as a variable. 

The question that arises is how robust the spin-interference effect is to deviations from the 
exact polygon-path. If we imagine a finite channel in the shape of a polygon in which the 
electron will undergo specular reflection on the walls, then in our model the path will be a 
zig-zag pattern and we might suspect that contributions to the rotation of the zig and zag 
paths will cancel each other as long as they are small. More specifically, if we choose the 
j/-axis as the direction of the rotation corresponding to motion along the polygon side and 
the direction of the zig and zag rotation under an angle ±/3 respectively with the y-axis in 
the x — y-plane, then a combined zig-zag gives a rotation 



where s' = V s 2 + w 2 with s the precession angle along the exact polygon side and with w 
the precession angle due to the finite width. If both the deviation angle (e$(0)) and the 
spin-orbit coupling (e^s')) are small then the correction due to the finite width is third- 
order. This applies also to a polygon with a large number of vertices and the deviation from 




(38) 



R e Rt = R z {(3)R y {s')R z {-2(3)R y {s')R z {l3) 




(39) 
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a circular trajectory will therefore be third order, explaining the rapid convergence of the 
interference pattern to that of a circle, as discussed above. 

If the spin-orbit coupling is finite then the finite width gives a linear correction to the 
rotation. Also, if the spin-orbit interaction is small and the deviation angle large, a linear 
correction results. Hence the calculated spin-interference effect is in general very sensitive 
to deviations from the exact polygon path and a finite width would result in an average over 
many possible paths that could potentially smear out the spin-interference. 

In summary, spin-orbit coupling in polygon-structures gives rise to spin-interference with 
strong dependence on the spin-orbit coupling constants and sidelength and alignment of the 
polygon. The dependence on the sidelength, alignment, and polygon- type can be exploited 
to extract the spin-orbit parameters from the spin-interference pattern which is observable 
for instance from its effect on the amplitude of the Aharonov-Bohm oscillation. 
The first of us would like to thank Harmen Warringa for reading the article. 
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FIG. 2: Interference intensity in the origin for the square as a function of the Rashba precession 
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angle sr and the linear Dresselhaus precession angle sdi in units of it. From bottom to top, the 
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Length (nm) Length (nm) 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Offset (tc/2) Offset (tc/2) 



FIG. 3: Interference intensity in the origin for the square as a function of the sidelength and 
alignment with the ic-axis for linear Dresselhaus coefficient am = 0.7 x 1CP 12 eVm and Rashba 
coefficient or = 1.5 x 10 -12 eVm for the left figure and a.R = 3.0 x 10~ 12 eVm for the right figure. 
L = 1000 nm corresponds to soi = 0.397T, sr = 0.847T respectively sr = 1.687T for an effective mass 
m* = 0.067 mo with mo the bare electron mass. The interval between contour lines is 0.25. 
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s s 
-4jc -2% R 2% 4% -4jc -2tt r 2k 4tc 



FIG. 4: Interference intensity in the origin for the triangle (upper left), square (upper right), 
hexagon (lower left) all aligned along the ai-axis and the circle (lower right). The sidelengths of the 
triangle and the hexagon and the radius of the circle have been scaled according to: 3L tr jangie = 
6Lhexagon = 2ttR = 4L square and the units of sr and sdi correspond to the rotation angle over one 
sidelength of the square. The interval between contour lines is 0.25. 
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